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ABSTRACT 

Several recent studies have shown that about half of the massive galaxies at z ~ 2 are in a quiescent 
phase. Moreover, these galaxies are commonly found to be ultra-compact with half-light radii of 
~ 1 kpc. We have obtained a ~ 29 hr spectrum of a typical quiescent, ultra-dense galaxy a,t z — 2.1865 
with the Gemini Near-Infrared Spectrograph. The spectrum exhibits a strong optical break and several 
absorption features, which have not previously been detected in z > 2 quiescent galaxies. Comparison 
of the spectral energy distribution with stellar population synthesis models implies a low star formation 
rate (SFR) of 1 - 3 M© yr'^, an age of 1.3 - 2.2 Gyr, and a stellar mass of - 2 x 10" Mq. We detect 
several faint emission lines, with emission-line ratios of [Nii]/Ha, [Sii]/Ha, and [Oii]/[Oili] typical 
of low-ionization nuclear emission-line region. Thus, neither the stellar continuum nor the nebular 
emission implies active star formation. The current SFR is < 1% of the past average SFR. If this 
galaxy is representative of compact quiescent galaxies beyond z = 2, it implies that quenching of star 
formation is extremely efficient and also indicates that low luminosity active galactic nuclei (AGNs) 
could be common in these objects. Nuclear emission is a potential concern for the size measurement. 
However, we show that the AGN contributes < 8% to the rest-frame optical emission. A possible 
post-starburst population may affect size measurements more strongly; although a 0.5 Gyr old stellar 
population can make up < 10% of the total stellar mass, it could account for up to ~ 40% of the 
optical light. Nevertheless, this spectrum shows that this compact galaxy is dominated by an evolved 
stellar population. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

The first massive, quiescent galaxies (> 10 "M(d) arose 
when the universe was on ly ~3 Gyr old (e.g.,lLabbe et al.l 
[2005t iKriek et all [2006 b: Willia ms et almOOl o r per- 
hap s even earlier (e.g., Bra mmer fc van DokkumI [20071: 
Mancini et al. 2009; Fontan a et al.ll2009l : lMobasher et al I 
2005t iWikhnd et al.n2008i) . Remarkably, these g alaxies 
already form a red sequence at z ~ 2.3 (Kric k et al.l 
l2008bf) . The relatively youn g ages 0.5 Gyr ) and po st- 
starburst spectral shapes (|kriek et al.ll2006R l2008bf ) of 
the z ~ 2.3 red-sequence galaxies suggest that a sig- 
nificant fraction of the stars have formed over a short 
timescale in an intense starburst. Sub- millimeter bright 
galaxies are possible candidates t o represent this vigor- 
ous phase of star formation (e.g., fChapman et al.l [20041 : 
[Tacconi et aL[[2008[: [Wall et al.[^008D . These dustv star- 
burst galaxies have observed star formation rates (SFRs) 
of several hundreds up to a thousand solar masses a 
year. The exact mechanism responsible for transforming 
such active syste ms into quiesce nt galaxies i s still sub- 
ject to debate (e.g..[Croton et al.l 2006: Bower et al.[[2006[: 
Naab et al. | [2007l: [Dekel &: Birnboimii2008(iHopkins et all 



escent systems at these e pochs fe.g.. [van Dokkum et al.l 
[20061 : [Tavlor et al1[2009af ). The majority of them quench 
or assemble into more massive systems at later times, 
and the number density of the massive end of the red 
sequence at z ^ 2 .3 is only ^ 1/8 of the local value 
([Kriek et al.[[2008bl ). Furthermore, the future evolution 
of massive quiescent galaxies at z ~ 2.3 is still un- 
clear. Their evolved stellar populations suggest that 
they passively evolve into their local analogs. How- 
ever, their strong size and slow color evolution contra- 
dict this picture. Recent morphological studies show 
that massive quiescent systems at z ~ 2 are remark- 
ably compact with effect i ve radii of ^ 1 kpc (e.g., 
Truiillo et aC '2006', '2007^, ' Zirm et al.[[2007t [Toft et al l 



2007t [Long hctti et al. 2007: [van Dokkum et al.[ \200A 
Cimatti et al.l i2 008: Franx et al. 2008; Damianov et alj 
20081: [van deTWel et al. .2008 : Saracco et al. 2009 ). Lo- 
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Not all local early types were already massive, qui- 
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cal early types o f similar stellar mass are about a factor 
of 5 larger (e.g.. [van Dokkum et al.[ [20081 ). Thus, these 
high-redshift galaxies must evolve significantly after z ~ 
2, probably by ins ide-out growth , primarily through mi- 
nor mergers (e.g.,.Bezanson et al.li2009:.Naab et al. 20091 : 
[van der Wei et al.l[2009l: [Hopkins et"al][2009b[ ). In addi- 
tion, the slow color evolution of the red sequence from 
z ~ 2.3 to the present implies that passive evolution 
alone cann ot explain the observed color-redshift relation 
([Kriek et al . 2008b). 

However, both the size and color evolution studies are 
hampered by many uncertainties and detailed, crucial 
information on the early phases of massive, early type 
galaxies is still lacking. The constraints on the stellar 
populations and SFRs at high redshifts are poor, even 
for galaxies with spectroscopy and mid-infrared photom- 
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Fig. 1.— NICMOS image of 1255-0. This imase is 3'.'8 b y 3'.'8 ; 
thus, 1" corresponds to 8.3 kpc. See Ivan Dokkum et al.l II2008I) 
for detailed information and for the images of a large sample of 
such compact quiescent galaxies at z ~ 2.3. The galaxy is barely 
resolved; most of the visible structure is due to the point-spread 
function. 



etry (e.g.. lKriek et al.ll2008aHMuzzin et al.ll2009| ). Thus, 
it is still unclear how "dead" z > 2 quiescent galaxies 
really are. Also, there are no dynamical mass measure- 
ments available for quiescent galaxies beyond z = 2. 
Consequently, all stellar mass estimates are photometric, 
and thus suffer from uncertainties in the derived stel- 
lar populations and from assumptions in the metallic- 
ity and the initial mass function (IMF). Moreover, spec- 
troscopic redshifts are extremely difficult to obtain for 
quiescent gala xies without emissio n lines. Using optical 
spectroscopy, Cimatti et al.l (|2008l ) have derived several 
absorption-line redshifts, but due to the relative faint- 
ness of quiescent galaxies in the rest-frame UV these ob- 
servations require ~ 100 hr of integration and result in 
incomplete samples. Near-infrared (NIR) spectroscopy 
allows the detection of the "bright" rest-frame optical 
continuum emission, but deep spectra are expensive due 
to the lack of multiplexing, the bright NIR background, 
and strong OH lines. Thus, so far no rest-frame opti- 
cal absorption lines have been detected in NIR spectra 
oi z > 2 quiescent galaxies, and previous redshift de- 
terminations rely on the detection of the Balmer and/or 
4000 A break with uncertainties of Az/(1 + z) < 0.019 
(iKriek et al.ll2008af) . 

In this paper we present a ^ 29 hr Gemini Near- 
Infrared Spectrograph (GNIRS: 'Eh as et al.ll2006[) spec- 
trum of a compact, quiescent galaxy at z = 2.1865, al- 
lowing a detailed study of the stellar population and the 
detection of any rest-frame optical emission and absorp- 
tion lines. Moreover, it gives us a glance into the fu- 
ture, as this is the deepest singlc-slit NIR spectrum ever 
taken of a 2 > 2 galaxy. Throughout the paper, we as- 
sume a ACDM cosmology with Q,^ = 0.3, JIa = 0.7, and 
Hq — 70 km Mpc~^. All broadband magnitudes are 
given in the Vega-based photometric system. 

2. TARGET SELECTION AND DATA 

The target was chosen from our GNIRS NIR spec- 
troscopic survey f or massive galaxies at z ~ 2.3 
(jKriek et al.l l2008af l. All galaxies were originally se- 
lected from the Multi- Wavelength Survey by Yale- 
Chile (MUSYC) , which provid es us with deep op tical- 
IR photometry (jGawiser et al.ll2006HQuadri et al.ll2007l : 
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Fig. 2. — Comparison of the broadband photometric rest- 
frame UV-NIR BED of 1255-0 (black open squar es) with the 
full galaxy sample (colored circles) from MUSYC IjQuadri et al.l 
l2007IMMarchesini etahl 120091) . To match the redshift and stellar 
mass of 1255-0, we restrict the sample to galaxies in the range 
1.9 < Zphot < 2.5 and 11.0 < log(M*/M0) < 11.5. The SEDs are 
normalized to unity at 1 /^m. The color of each SED reflects the 
rest-frame UV flux in the normalized spectrum. Triangles indicate 
1(T upper limits. The gray curve represents the best-fit BPS model 
to the full BED of 1255-0. 



Tavlor et al]l2009bHDMnen et al.l[200l iMarchesini et all 
20091 ). These "shallow" spectra (typically ~3 hrs of in- 



tegration, see iKriek et"a l. 2008a) allowed us to derive 
continuum redshifts and classify the galaxies. 

We selected 1255- from the nine ma ssive quiescent 
galaxies presented in IKriek et al.l (|2006bf) . because of its 
redshift (the optical break falls in the J band), and its 
visibility at the time of our GNIRS run. The galaxy is 
not brighter than the other candidates. Also its effective 
radius (rg) of 0.78 kpc is very similar to the median 
(0.9 kpc) o f the other massive, quie scent galaxies in our 
sample (see van Dokku m et al.l2008[ ). Figure[T]shows the 
image of 1255-0 as obtained by the Near-Infrarcd Camera 
and Multi-Object Spectrometer (NICMOS). Altogether, 
this galaxy seems typical for the general class of quiescent 
galaxies at z ^ 2.3. 

We also consider how this galaxy compares to the gen- 
eral population of massive galaxies at similar redshift 
(including star- forming galaxies). In Figure [5] we com- 
pare the rest-frame UV-NIR spectral energy distribu- 
tions (SEDs) of all galaxies in the deep MUSYC sample 
(|Quadri et all l2007t IMarchesini et all I2OO90 with com- 
parable redshift and stellar mass. The redshifts are all 
photo metric, and derived using EAZY (H rammer et al.l 
l2008f l. The stellar masses are d erived using the code 
described in the Appendix, for the lMarastonI ()2005l ) stel- 
lar population syn t hesis (SPS) models, solar me tallicity, 
the iCalzetti et ahl (j2000f ) reddening law, and a iKroupal 
12001') IMF. Figure H] illustrates that 1255-0 is slightly 
redder than the average massive galaxy at this redshift. 
This is expected given the fact that it was selected to be 
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TABLE 1 

Observations 







Exposure Time 


Seeing 


ID run 


Date 


(minutes) 


(") 


GS-2005A-Q-20 


2005 May 19 


85 


0.9 




2005 May 27 


80 


1.0 




2005 May 30 


80 


1.0 


GS-2006A-C-6 


2006 Feb 24 


60 


0.5 


GS-2007A-C-9 


2007 Mar 11 


360 


0.5 




2007 Mar 12 


330 


0.4 




2007 Mar 13 


340 


0.5 




2007 Mar 14 


390 


0.6 



quiescent. 

In total we integrated nearly 29 hr on 1255-0 with 
GNIRS (with individual exposures of 10 minutes) divided 
over three observing runs in 2005 May, 2006 February, 
and 2007 March. The integration times and average see- 
ing are given in Table [T] During the first three nights the 
conditions were mediocre, with occasional clouds and an 
average seeing of 1". The weather conditions were ex- 
cellent during the last two runs, with clear skies and an 
average seeing of 0.5". 

The galaxy was observed in a cross-dispersed mode, 
in combination with the 32 lines mm~^ grating and the 
0'.'675 slit. The spectral resolution varies between ^--^900 
and ~1050 over the different orders. Observations were 
done using an ABA'B' on-source dither pattern, such 
that we can use the average of the previous and fol- 
lowing exposures a s sky frame (jvan Dokkum et"al]l2004l : 
iKriek et al.l l2006al ). Acquisition was done using blind 
offsets from nearby stars. Before and after every observ- 
ing sequence we observe an AVO star, for the purpose 
of correcting for telluric absorption. The final spectra of 
the two stars are combined to match the airmass of the 
observing sequence. 

A detailed description of the reduction procedure of 
GNIRS cross-dispersed spectra is given in iKriek etldl 
(|2006aD . In summary, we subtract the sky, mask cosmic 
rays and bad pixels, straighten the spectra, combine the 
individual exposures, stitch the orders, and finally cor- 
rect for the response function. The different observing 
sequences are weighted according to their signal-to-noise 
ratio (S/N) when being combined. 

A one-dimensional spectrum is extracted by summing 
all adjacent lines (along the spatial direction) with a 
mean flux greater than 0.1 times the flux in the central 
row, using optimal weighting with the S/N. We extract 
both a high- and low- resolution spectrum with 10 A and 
50 A per bin, respectively. The high-resolution spectrum, 
which is resampled such that no resolution is lost, is used 
for spectral features, while the low-resolution spectrum 
is extracted to study the continuum emission. In order to 
flux calibrate the spectrum, we derive the spectroscopic 
J, , and K fluxes by integrating over the corresponding 
filter curves. We derive one scaling factor by comparing 
the spectroscopic colors with the J, , and K broadband 
photometric data from MUSYC. Finally, we multiply the 
spectrum by this scaling factor. 

The low-resolution spectrum is shown in Figure [3l 
Note the clear detection of the rest-frame optical break 
in the J band and the relatively high S/N in the con- 
tinuum (compared to the spectrum shown in IKriek et al.l 



I2006bf) . 

3. ANALYSIS 

We study the properties of this galaxy in two ways, first 
by measuring and analyzing the spectral features (Sec- 
tions 13.11 and 13. 2|) and second by modeling and decom- 
posing the stellar continuum emission (Sections 13. 3I - IX5)) . 
Finally, in Section 13.61 we compare the deep spectrum 
modeling results with our previously published shallow 
spectrum. 

3.1. Emission and Absorption Features 

We measure the redshift and all emission-line prop- 
erties by modeling the extracted high-resolution (10 A 
per bin) one-dimensional spectrum. We detect Ha, [Nil] 
AA6583, and [Sii] AA6717 in the K band. The bottom 
panels of Figures [4] and [5] show the relevant part of the 
two-dimensional and one-dimensional spectrum, respec- 
tively. For the former we first removed the best-fit con- 
tinuum model (see Section 13. 3p to make the lines more 
visible. We detected none of these lines in our shallow 
spectrum (jKriek et al.ll2006bh . as they are too faint. 

We model Ha and the two [Nil] lines simultaneously, 
by assuming the same redshift and width for all three 
lines, the best-fit continuum model as derived in Sec- 
tion l3.3[ and Gaussian profiles. As the detected emission 
lines are all faint, it is important to accurately correct 
for continuum emission (including the Balmer absorp- 
tion lines). The continuum model is corrected for the 
spectral resolution of GNIRS and convolved to the same 
velocity width as the emission-line model. Furthermore, 
we adopt the ratios of transition probability between the 
two [Nil] lines of 0.34. 

We derive the best values for the redshift, line width, 
and the fluxes of the emission lines by minimizing . 
The uncertainties on the modeling results are derived 
using 500 Monte Carlo simulations. In the simulations, 
we perturb the measured spectrum using the noise spec- 
trum. The results are listed in Table Si [S ii] AA6717 
and [Sil] AA6731 are modeled in the same way; thus the 
redshift, width, and scaling for both lines were free pa- 
rameters. 

In the H band we detect no emission lines (see Figure 
[5] top-right panel). We derive 2 a upper hmits using 500 
Monte Carlo simulations in which we model the lines 
assuming the width and redshift as obtained from the 
Ha and [N ii] lines, and the continuum model. We fit all 
three expected lines H/3, [O ill] AA4959 and [O ill] AA4959 
simultaneously, assuming a ratio of 0.33 between the two 
[Oiii] lines. Upper limits are derived from the best-fit 
results of all simulations. 

In the J band we detect the doublet [O ii] (see Figures 
[3]and[51 top- left panels). We do not resolve the two lines 
separately, and thus we fix the ratio (to 1) and redshift 
when fitting the lines. The combined fitting parameters 
are given in Tabled 

For each pair of lines we fit simultaneously, we assume 
that all line emission originates from the same physi- 
cal region. [Sii] and Ha yielded consistent redshifts, 
but the different combinations of lines resulted in dif- 
ferent widths. This might imply that the lines originate 
from different processes and regions in the galaxy. How- 
ever, the lines are all faint, and the errors on the line 
widths are probably somewhat underestimated because 
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Fig. 3. — The ~ 29 hr deep GNIRS spectrum of 1255-0 (black). The spectrum is sampled to bins of 50 A in observed frame, and 
smoothed by a boxcar of 75 A. The gray shaded background indicates the noise leve l of the spectrum, w i th da rk being noisier. We omitted 
the regions in between the atmospheric windows. Overplotted in green is the best- fit [Bruzual fc CharlotI l|2003l ) SPS model to the spectrum 
and optical-IR photometry, assuming a lChabrieil 1)20031 1 IMF and solar metallicity. This fit corresponds to a stellar mass of 2.3 X 10" M0, 
an age of 2.1 Gyr, a t of 0.3 Gyr, a reddening of = 0.25 mag, an SFR of 1.9 Mq/yv, and a specific SFR of 0.008 Gyr"!. The locations 
of possible absorption and emission lines are indicated by the red dashed and blue dotted lines, respectively. 



of our fitting procedure (allowing more freedom in the 
fits would result in larger uncertainties, particularly in 
the linewidths). Furthermore, we do not know the rel- 
ative contributions from the two [O ii] line. Thus, we 
cannot draw any firm conclusion from the different line 
widths. 

As far as we are aware of, this is the first rest-frame 
optical spectrum of a quiescent galaxy beyond z = 2 for 
which rest-frame optical absorption lines are detected. 
Figure [6] shows that we detect Hry, H^, Caii-K and Caii- 
H/He. Caii-K and Caii-H/He are also visible in the 
top-right p anel of Figure IH In Figur e [6l we also show 
the best-fit iBruzual fc CharlotI ()2003l ) models for fixed 
solar metallicity and when leaving the metallicity as a 
free parameter (see Section . Furthermore, the ab- 
sorption lines Mg b A5175 and H/3 are detected in the 
low-resolution spectrum in Figure [3l There might also 
be a hint for H(5 in Figure [3l However, a significant de- 
tection of this line requires an even deeper spectrum. 

This deep spectrum allows an accurate measure of the 
strength of the 4000 A break, which is an indicator of 
the evolution stage or age of the s tellar population. We 
use the definition of iBalogh et al.l (119991 ) to measure the 
strength of this break, and find a value of -D„(4000) = 



formation history (SFH), the metallicity, and the dust 
content of the galaxy. For the extreme case of a simple 
stellar population (SSP) model with solar metallicity, an 
age of ^ 0.6 Gyr is needed to produce such a break. 
Assuming little or no dust, this value can be seen as 
a lower limit on the population age. A more detailed 
measurement of the age of the stellar population will be 
given in Section [ 



1.40 



+0.03 
-0.03- 



The corresponding age depends on the star 



3.2. Ionization Mechanism 

Emission-line ratios can be used to study the origin 
of the ionized emission. In particular, [Nil] /Ha, which 
reflects both the metallicity and ionization parameter 
C/ of a galaxy, is a powerful discriminator. Ratios of 
[Nil] /Ha > 1 suggest that an object is ionized by a hard 
radiation fie ld; Hli regions are not able to produce such 
high ratios (jKewlev et al.ll2001h . The ratio of [Sii]/Ha*^ 
can be used for a similar purpose, as Hii regions have 
[Sil]/Ha < 0.6. The high values for both [Nii]/Ha and 
[Sii]/Ha of l-itl'l and 0.8j:o;2 respectively show that 
normal star-forming regions are not the dominant con- 
tributor to the line emission in this galaxy. 

^ For [S ll] we take the sum of both lines, while for the ratio of 
[Nil] /Ha only [Nil] AA6583 is used 
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TABLE 2 

Spectral Features 





A 




EWs 


F 


L 




ad 


Lines 


A 




A 


10~i^ ergs s~i cm^2 


10''^ ergs s~i 


z 


km s^i 


[O ll]^ 


3727 




14.81^;^ 


9 9^+0.31 

^•^^-0.36 


n oT+O-ll 
0-87Io.i3 






M/J • ' 


4ooo 




< i.O 


< U.4Z 


< U.IO 






\t~\ tttI a.b 
[\J IIIJ ' 


4yoy 




< U.4 


^ n 1 y1 
< U.14 


< U.Uo 






[O III]''''' 


5007 




< 1.3 


< 0.43 


< 0.15 






[Nil] 

[Nil] 
[Sii] 
[Sii] 


6548 
6563 
6583 
6717 
6731 




1 7+0.1 
^■'-0.1 

4 0+0-5 

5 1+0--* 

9 O+0.2 

^"-"-o.a 
9+0-3 


55+0-05 
'^-''''-0.03 

1 oq+o-w 

^-^^-0.18 


20+° 02 
"^■^"^-0.01 
46+0 06 

60+0-05 
'-'■"'-'-0.03 
n 97+0.03 
"•^'-0.03 

10+0-04 
'-'■l''-0.03 


9 lQfic;+0.0002 
9 iofic:+0.0002 
9 lQfiK+0.0002 
9 1OR4+0.0003 
9 1OR/I+0.0003 


2031?^ 
203li^ 
203t?^ 

108lg 


Fixed to 


best 2spcc as 


derived 


from Ha and [Nil] 


1* Fixed to best cr as derived 


from Hq and [Nil] 


Corrected for Balmer 


absorption.^ The 


resolution c 


f GNIRS is R 


~ 1000, 


corresponding to a 


= 128 kms^i 









TABLE 3 
Stellar population properties 





Log M, 


Log age 


Log T 


Av 




Log SFR 


Log SFR / M, 


UBVRIz + NIR spectrum 


(Mq) 


(yi) 


(yi-) 


(mag) 


z 


(Afeyr-i) 


(Gyr-i) 


BC03, Chabrier, Zq 
BC03, Chabrier, Z free 
BC03, Chabrier, Zq, +IRAC 
BC03, Chabrier, Z free +IRAC 
MA05, Kroupa, Zq 
MA05, Kroupa, Zq +IRAC 


11 26+°°° 

il.ZD_Q Q2 

11 36+°°° 
11 26+°°° 

il.ZD_g 02 

, , 9-+0.00 

, 94+0.00 
11-^^-0.04 


Q 34+0.02 
^-'-'^-0.04 
q 10+0-02 
y-l'-'-0.02 
Q 09+0.04 
^-'-'^-0.04 
q 10+0-02 
y-l'-'-0.02 

q 26+°-°° 
^-^"-0.02 

q 26+°-°° 
^-^"-0.04 


8.51°:? 

8.31°:? 
8.51°:? 
8.31°:? 

8 4+0-0 
^-^-0.1 

8 4+0-0 

^-^-0.1 


20+°-l5 

30+° °^ 
'-'-■''^-0.05 

25+0-10 

'J-^^-o.io 

30+° °° 
'-'-■''^-0.10 

00+° °5 

00+°-^° 
i-'J'^-o.oo 


02+° °° 
'J-'J^-o.oo 

05+°-°° 

02+° °° 
'-'-'J^-0.00 

05+°-°° 

02+° °° 
'J-'J^-o.oo 

02+° °° 
i-'J^-o.oo 


19+°-l<3 
"•l='-0.20 

o.45lg:°° 

0.281°:°^ 
45+°-°° 

OQl + 0.11 
'J-'-'l-0.42 

QQ-,+0.20 

'J-'-'l-0.42 


9 17+0.17 
^■l'-0.18 
1 O1+0.00 

-l-»l-0.34 
9 n«+0-08 

-^.U»_Q 27 
1 O1+0.00 

-1-O1-0.34 
9 90+0.12 
^-^•^-0.40 
9 90+0.21 
^-^•^-0.40 



Note. — Wlien deriving the confidence intervals we only use the model grid points enclosed by the 68% x threshold values (see Appendix) 
Due to this discretization, non-zero confidence intervals may be rounded to zero. 



The ratio of [O iii]/[0 ii] can be used to further charac- 
terize the hardness of the ioni zing radiation field (jShieldsl 
\Tmi iKewlev fc DoDital[200l . [Oiii] AA5007 is not de- 
tected in our spectrum, but we derive a, 2 a upper limit on 
this parameter. We find log([Oiii]/[Oiip < -0.89, im- 
plyin g an ionization parameter U < 10* (|Kewlev et alj 
|2001|) . This ratio is consistent with the original definition 
for low-ionization nuclear emission-line regions (LINERs) 
of [Oii]/[Oiii]> 1, and so the line emission in 1255-0 is 
most likely caused by a LINER. 

The origin of LINER emission has been a subject of de- 
bate since its original classification as a n active galactic 
nucleus (AGN) class bv lHeckman l (1980). Although there 
is substantial evidence that many LINERs are powered 
by accretion onto massive black holes, LINER emission 
can also originate from a young starburst or by shock 
heating through cloud collisions in duced by galaxy merg- 
ers or starburst-driven winds (e.g.. lDopita fc Sutherland! 
[1995). Because the SFR in 1255-0 is very low (see Sec- 
tion [3?3l), a starburst or starburst-driven wind is unlikely 
to cause the LINER emission, and thus a low-luminosity 
AGN is the more likely option. Nonetheless, even if the 
observed LINER emission does not originate from black 
hole accretion, it is still the case that normal star-forming 
regions do not dominate the observed line emission. 

Finally, we note that the line widths are not necessarily 
indicative of the depth of the gravitational potential, as 
in the local universe there is a large scatter between the 
gas line widths of LINERs and the velocity dispersi on of 
the stars of their host galaxies (| Greene fclfol[2005h . 



3.3. Modeling the Continuum Emission 

We study the nature of the stellar population by 
comparing the sp e ctrum wit h the SPS mo dels of 
iBruzual fc Chariot! (|2003f ) and llVIarastonI (|2005f) . using 
the code described in the Appendix. We assu me an ex- 
]X)nentially declinin g SFH with timescale Tj_ a Chabrie 
2003) or Kroupal (|200H ) IIVIF, and the iCalz ctti ct al 
2000) reddening law. We derive model spectra for a grid 



with T between 10 IVlyr and 10 Gyr in steps of 0.1 dex, 
Ay between and 3 mag in steps of 0.05 mag, and age in 
steps of 0.02 dex with a minimum age of 10 IVIyr and the 
maximum age not exceeding the age o f the universe. The 
metallicity when fitting the [Maraston (2005) model is 
fixed t o solar {Z = 0.02), but for the IBruzual fc CharlotI 
(|2003f) models we vary the metallicity, among subsolar 
{Z = 0.004), solar, and supersolar (Z = 0.05). The red- 
shift is fixed to z = 2.1865 as derived from the emission 
lines (see Section 13.11) . We fit the low-resolution one- 
dimensional spectrum (50 A per bin in observed frame) 

and search for the best solution by minimizing y^. 

In contrast to our previous studies (Kriek et al.l 
I2006airb . 2008a ). we do not mask regions with low atmo- 
spheric transmission or strong sky emission when fitting 
the spectrum. As the S/N of this spe ctrum is consider- 
ably higher than in lKriek et al.l (j2008af ). we apply a much 
smaller bin size. Bins in bad wavelength regions will have 
larger uncertainties, and thus simply have l ower weight 
in the fit. When using larger bins (400 A in iKriek et al.l 
l2008a^ this method is less appropriate, as a bad region 
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Fig. 4. — Selected regions of the two-dimensional spectrum. 
Each pixel is 25 A in observed frame (8 A in rest-frame) in the 
wavelength (horizontal) direction by 0'.'15 in the spatial (vertical) 
direction. The top-left, bottom-left, and bottom-right panels show 
the emission lines [O ll] , Ha and [N ll] , and [S ll] , respectively. The 
continuum emission has been removed, using the best-fit model, to 
make the lines more visible. The top-right panel shows the wave- 
length region around the break. Call-H and Call-K are clearly 
visible in the two-dimensional spectrum. In this panel, the contin- 
uum emission has not been removed. 



will contaminate nearly all bins. 

To further constrain the SFR, we include the rest- 
frame UV broadband photometry. Furthermore, we ex- 
tend the SEP into the rest -frame NIR using IRAC pho- 
tometry ()Marchesini et al.l 120091. The rest-frame NIR 
helps to constrain the stellar mass and age of the galax;^^ 
(e.g. j Labbc et al. 2005 



[200 



Shapley et al. 2005; Wuvts et alj 
Muzzin et al.ll2009D . However, as the SPS models 
are still uncertain in this regime, we fit the galaxy both 
with and without the IRAC photometry. 

We derive 68% confidence intervals on all stellar pop- 
ulation properties using 200 Monte Carlo simulations, 
as described in the Appendix. The photometric un- 
certainties are increased using the template error func- 
tion, which accounts for uncertainties in the model 
templates as a function of the rest-frame wavelength 
(p3rammcr ct al. 2008). Furthermore, we apply the au- 
tomatic scale option, such that for each simulation the 
simulated spectrum is calibrated using the simulated J, 
H, and K photometry (see Section[2]and the Appendix). 
Thus, the uncertainty in the calibration of the spectrum 
is explicitly taken into account. 

In Table [3] we give all modeling results for the dif- 
ferent SPS libraries, free or fixed metallicity, and with 
or without IRAC. In Figure [3 we s how the spectrum 
and best-fit iBruzual fc Chariog (|2003D model with solar 
metallicity. The continuum fitting implies a stellar mass 
of ~ 2 X 10" Mq, a reddening of Ay = 0.0 - 0.3 mag, 
a star formation timescale of r ~ 0.3 Gyr, an age (since 
the onset of star formation) of 1.3 — 2.2 Gyr, and an SFR 
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Fig. 5. — One-dimensional spectrum (10 A per bin in ob- 
served frame) in regions around the emission lines. Fx is in 10^^^ 
ergs~^ cm~^ A~^. The red and blue lines represent the best con- 
tinuum and emission-line fits, respectively. While we detect [On], 
Ha, [Nil] and [Sll], we see no signs of H/3 in emission or the two 
[Olll] lines. 
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Fig. 6. — Spectrum (20 A per bin in observed frame) in the 
wavelength region around the rest-frame optical break. Fx is in 
10"^® ergs^^ cm~-^ A~^. The red and blue curves represent the 
best fits for supersolar (2.5Zq) and solar metallicity, respectively. 
Both the Call-H/He and Call-K lines are detected, as well as Hjj 
and HC. 



of 1 — 3 Mq yr ^. The universe is 3 Gyr at z ~ 2.1865, 
which implies a formation redshift Zform = 4 — 7. 
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Fig. 7. — Reduced vs. continuum redshift (number of degrees 
of freedom N^^f = 31 9), when leaving; redshift as a free parame- 
ter in the SPS model IjBruzual fc Charlotj[2003 ) fitting. The other 
parameters, including metallicity, are allowed to vary. The rest- 
frame UV-NIR broadband photometry is not included in the fit. 
The 68% and 95% confidence intervals are indicated by the dotted 
lines. The best-fit emission line redshift is indicated by the ver- 
tical solid line, and its 68% and 95% confidence intervals by the 
green and blue shaded regions, respectively. The continuum emis- 
sion yields a well-constrained spectroscopic redshift. The redshift 
difference between the stellar and nebular emission is not at all 
significant. 



It is remarkable how well constrained the formal con- 
fidence intervals are. However, the formal errors do not 
reflect the true uncertainties properly, as they are domi- 
nated by the systematic effects such as the assumptions 
concerning the SPS mo dels, the SFH, metal l icity, and ex- 
tincti on law (see e.g.. IShaplev et aL 20 011: IWuvts et al] 
20071: iKannappan fc Gawiseil l2007t iConrov et all 120081 
Muzzin et al.l l2009l for mor e discu ssion on this topic). 



For the iBruzual fc Charlotl ([2003) models, supersolar 
metallicity is formally preferred above solar. We do not 
believe the fact that this result is significant, due the 
strong degeneracy between age and metallicity, and the 
uncertainties in the SPS model, and so the metallicity 
could well be less, e.g., solar. Nonetheless, the different 
models are surprisingly consistent. The inclusion or ex- 
clusion of IRAC has little effect on the m odeling results 
for this deep spectrum (see lMuzzin et al.|[2009i . for a gen- 
eral discussion of the inclusion of IRAC data in the mod- 
eling of z ~ 2.3 galaxies with shallower spectroscopy). 

The detection of the strong break and absorption fea- 
tures allows the measurement of a stellar continuum red- 
shift. We use the same fitting procedure as described 
above; only this time, we leave the redshift as a free 
parameter. Furthermore, we fit only the spectrum; thus, 
the rest-frame UV and NIR broadband photometric data 
are not included. In Figure [3 we show the reduced 
as a function of redshift. We find a stellar redshift of 
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The emission-line redshift is also indi- 



observed 

Fig. 8. — Multi-component modeling of 1255-0. The black solid 
dots and blue SEDs show the original broadband photometric data 
and best-fit SPS model to the full SED (including NIR spectrum). 
In the top panel we assess the contribution from a post-starburst 
population (with an Ay of 0.25 mag), which can account for < 10% 
of the stellar mass and < 40% of the i?-band flux. The corrected 
SED and corresponding best-fit SPS model are represented by the 
open symbols and white curve, respectively. The inset shows the 
quality of the fit as a function of the contribution of the post- 
starburst population to the total stellar mass. The red lines indi- 
cate the maximum contribution. In the bottom panel, we assess a 
contribution from an AGN. The inset shows the quality of the fit 
as a function of the AGN contribution to the H-hand flux. The 
green lines indicate the maximum contribution (8%). In both cases, 
these flts imply that the galaxy is primarily composed of an evolved 
stellar population. 

cated in this figure. There is no evidence for a significant 
redshift offset between the stellar and nebular emission. 
Nonetheless, we cannot exclude a possible offset, due to 
the relatively large uncertainties on the continuum red- 
shift. 

3.4. A Recent Starburst? 

When modeling the SED by SPS models, we assume an 
exponentially declining SFH. However, this is a simplifi- 
cation, and more complex SFHs are more realistic. For 
example, subsequent merging is expected to result in cen- 
tral dissipational starbursts (e.g., Hopkins et al. 2009a). 
This would result in a younger stellar population in the 
central part of the galaxy, with a lower mass-to- light ratio 
(M/L) than the older underlying population. Similarly, 
if the galaxy experienced recent star formation in a disk- 
like component, the galaxy is composed of an old central 
concentration and a more extended young component. 
This younger population will contribute relatively more 
to the observed light than to the stellar mass. 



8 



Kriek et al. 



In order to assess whether a recent starburst took place 
and how much it may contribute to the stellar mass of 
the galaxy, we investigate how much of the light can be 
accounted for by a 0.5 Gyr SSP model, with 0.25 mag 
of visual extinction (see Table [3]). We choose a post- 
starburst instead of an ongoing starburst, as we know 
that the current global SFR in the galaxy is very low. 
We subtract different mass contributions and apply the 
same fitting procedure as discussed in Section [5T51 The 
maximum contribution of the post-starburst population 
is set by the rest-frame UV flux, and shown by the red 
SED in the top panel of Figure El Up to -10% of the 
stellar mass may have been formed in a recent starburst. 
Due to the relatively low M / L of this post-starburst pop- 
ulation, the contribution to the H-h&nA light is much 
larger 40%). 

In the inset of the top panel of Figure [5] we show the 
reduced of the best-fit SPS model to the corrected 
spectrum and photometry, for different mass fractions 
of the post-starburst stellar population. Models with a 
mass fraction between ~ and 5% provide equally good 
fits, while higher mass fractions result in worse agree- 
ment. The red line indicates the maximum contribution, 
corresponding to the red SED. Thus, while this galaxy 
may have experienced a recent starburst, both the light 
and the stellar mass are dominated by an older stellar 
population. Nonetheless, it still remains to be explored 
how much the light distribution may be different from 
the stellar mass distribution in these compact quiescent 
galaxies. 

3.5. Continuum Emission from an AGN? 

In Sect ion [3?2l we found that the line emission of 1255-0 
is of LINER origin. This raises the question of whether 
an AGN may contribute to the continuum emission of 
1255-0. We investigate this by subtracting different AGN 
contributions and fitting the corrected SED by SPS mod- 
els. We assume a power-law SED for the AGN. The 
maximum AGN contribution is set by the rest-frame UV 
fluxes in combination with the 8 /im IRAC band, and 
shown by the green line in the bottom panel of Figure [H 
The corresponding contribution to the H-hand flux is 
-8%. 

In the inset of the bottom panel of Figure [8] we show 
the reduced of the best-flt SPS model for different 
assumed AGN contributions. The green line indicates 
the maximum contribution, corresponding to the green 
SED. The flt clearly worsens for an increasing AGN con- 
tribution. Altogether, an AGN is unlikely to contribute 
signiflcantly to the continuum emission, and is limited to 
a maximum of — 8% to the H band. 

3.6. Comparison to the shallow spectrum 

While our 5 h r^ shafl ow spectrum of 1255-0 presented 
in iKriek et al.l (|2008bl ) provided a similar stellar mass 
and SFR (when corrected for the difference in the as- 
sumed IMF), the age, the dust content, and the con- 
tinuum redshift differ by ~ 2 cr compared to the 29 hrs 
deep spectrum. Our shallow spectrum of 1255-0 yields 
a continuum redshift of 2.3llg'Q7, an age of 0.57^Q2g 
Gyr, and an Ay of 1-2^0 e mag. Thus, this galaxy was 

The weather conditions during the first runs were significantly 
worse, and the effective exposure time is closer to 1-2 hr. 



previously classifled as a dusty post-starburst galaxy, at 
slightly higher redshift. In the fit to the shallow spec- 
trum the dominant optical break was thought to be the 
Balmer jump, while in our deeper spectrum the 4000 A 
break is found to be the more prominent one. If we fit 
the shallow spectrum of 1255-0 with the redshift fixed to 
z = 2.1865, we find a significantly older age of 2.9 Gyr 
and a lower Ay of 0.35 mag. 

The change in redshi ft also influences th e rest-frame 
color determination fsee lKriek et al.l l2008bV 1255-0 was 
among our reddest galaxies, with a rest-frame U—B color 
of 0.36lg;Qg mag. The deep spectrum yields an U — B 
color of 0.52 mag . However, if we use the same method 
as in iKriek et al.l ([2008b) , thus determining the color of 
the best fit, we find U — B = 0.30. This difference is 
likely caused by the discrepancy between the spectrum 
and the flt around 1.15 /.tm. 

We do not expect the full sample of shallow spectra 
to suffer as severely from this degeneracy b etween red- 
shift and stellar population as 1255-0. In ' Kriek et al.l 
(2008a) we apply the same continuum fltting procedure 
to the emission-line galaxies in our sample, and flnd a 
good agreement between the emission-line and contin- 
uum redshifts. This sample contained several galaxies 
with SED shapes similar to those without emission lines. 
Nonetheless, while it seems more likely that 1255-0 is 
the largest outlier 2a), caution is required. This case 
illustrates that more deeper spectra are needed. 

4. THE STAR FORMATION ACTIVITY IN 1255-0 

Comparison of the stellar continuum emission with 
SPS models (Section l3.3p conflrms our earlier results that 
the SF R in 1255-0 is strongly suppressed (jKriek et al.l 
l2006bD . Depending on the SPS hbrary and the assumed 
metallicity, the best-fit SFR is 1 — 3 M0yr~^. As ex- 
pected for a galaxy with low-level star formation, the 
total dust extinction is low with values of Ay of 0.0-0.3 
mag. 

The low SFR is independently confirmed by the 
emission-line diagnostics. In the previous section, we 
noted that star-forming regions cannot be the dominant 
contributor to the line emission. Nevertheless, if we as- 
sume that the detected Ha emission is caused by just 
star formation, we can use the calibration 

SFR (Mgyr-i) = 7.9 x lO'^^ Lho (ergs"') (1) 

as given by Kcnnicutt (1998) to derive the SFR from 
the Ha luminosity. Assuming no dust extinction, 
the observed luminos ity woul d resu lt in an SFR of 
A.OtoiMQjT-^ for a iSalpeted (fl955l) IMF. This value 
w ould decrease b y a factor of ~ 1.8 for a lChabrieil (j2003l ) 
or lKroupal (|2001h IMF, but increases by a factor of 2.3 
if there is 1 mag of visual extinction in the star forming 
regions. Although we have a good constraint on the to- 
tal attenuation of the galaxy from the spectral modeling 
{Ay = 0.0 — 0.3 mag), it is difficult to estimate the ex- 
tinction in the line-emitting regions without measuring a 
Balmer decrement. The combination of Ha and the 2 a 
upper limit on H/3 gives a 2 cr lower limit on the ratio of 
Ha/H/3 of 3.07. This limit is close to the intrinsic ratio 
for Hii regions of 2.76, and thus sets no constraints on 
the dust content. Nonetheless, even when assuming that 
the line emission originates just from star formation, and 
assuming 1 mag of dust extinction in the H ii regions, the 
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obtained SFR is still only ~ 4 Mq yr ^ (for a Chabrier 
IMF). 
The 

This value is slightly dependent on the choice for the SPS 
model, and on whether metallicity is assumed to be solar 
or left as a free parameter. If left as a free parameter, 
supersolar is preferred over solar, although we argue in 
Section [3731 that this result is not significant. The best-fit 
star formation timescale is 0.2 — 0.3 Gyr, and the best- 
fit stellar age is 1.8 — 2.2 Gyr and ^ 1.3 Gyr for solar 
metallicity and supersolar metallicity, respectively. The 
age of the universe at this redshift is 3.0 Gyr and the 
formation redshift is Zform ~ 4 — 7. 

We study the evolution stage of the galaxy by means of 
the quenching factor (/sf , which measures by what factor 
the SFR has been reduced compared to the past average 
SFR. We define Qsf as follows: 



SFRcu 



SFR, 



= 1 



SFR, 



past 



M, 
age 



(2) 



In this equation is not the current stellar mass, but 
all the mass formed in the galaxy. This quenching factor 
is related to the ratio of the age to the star formation 
timescale t (i.e., how many e- folding times have been 
passed since the galaxy was formed) in the following way 



gsf = 1 



age 



1 - e- 



(3) 



Because age/r is a relatively well-constrained parame- 
ter, Qsi is a fairly robust m easure for the evolution stage 
of a galaxy. We note that iDave) (|2008f ) defined a com- 
parable factor: the star formation activity parameter 
agf = (M*/SFR)/(i// — IGyr). Instead of using the best- 
fit age of the galaxy, this parameter uses the Hubble time 
(tn) minus 1 Gyr. 

1255-0 has a q^f of Qmito^^^Ji- Thus, the SFR in this 
galaxy has been reduced by more than 99% and the 
current SFR is less than 1% of the average past SFR. 
This implies that the star formation has been strongly 
quenched since the major star formation epoch. Fur- 
thermore, in Section 13.41 we found that a 500 Myr SSP 
accounts for <10% of the stellar mass, which implies that 
at least 90% of the stellar mass is in stars older than 0.5 
Gyr, and thus formed beyond z = 2.6. 

5. DISCUSSION AND CONCLUSIONS 

Due to the introduction of large photometric sur- 
veys with deep N IR imaging (e.g., iLabbe et all l2003l : 
iFranx et al.|[200l . many quiescent massive galaxies be- 
yond z = 2 have been identified i n the past few yea rs 
(e.g..lForster Schreiber et al.ll2004l : lDaddi et al.ll2004allbl : 
ILabbe et al.ll2005n . Moreover, these galaxies are typically 
found to be ultra-compact, with stellar densities that are 
about 2 orders of magnitude large r than in local early 
type galaxies of similar mass (e.g iTruiillo et al. Il2006l 



mm 



20071 : 



Zirm et al.ll2007t iT oft et aLH2007HLonghetti et ail 
van Dokkum et al . 2008; Cimatti et al.ll2008l ) 



Follow-up spectroscopic studies have tried to verify 
the broadband photometric redshifts and stellar popu- 
lation properties of these ultra-dense quiescent galaxies. 
This turned out to be extremely difficult, and with opti- 
cal spectroscopy tens to hundreds of hours are required 



(jPaddi et al.ll2005l : ICimatti et al.ll2008[ ). due to their rel- 
ative faintness in the rest-frame UV. Evolved galaxies 
beyond z = 2 are much brighter at NIR wavele ngths, 
corresp onding to the rest-frame optical. In Krick et al.l 
flOOeU), we used the Balmer and / or 4000 A breaks 
to obtain redshift estimates for a sample of nine massive 
quiescent galaxies. However, exact redshift measurement 
from rest-frame optical absorption lines remained out of 
reach until this paper. 

By integrating for nearly 30 hr with GNIRS we suc- 
ceeded in detecting for the first time rest-frame optical 
absorption lines in an NIR spectrum of a compact, qui- 
escent galaxy beyond z = 2. This deep spectrum has 
full NIR coverage (^ 1.0-2.4 /zm). In addition to the ab- 
sorption features H?7, HC, Caii-K, Caii-H/He, H/3, and 
Mg b, we have detected [On], Ha, [Nil], and [Six] in 
emission. All emission lines are faint with luminosities 

of 0.1 0.9 X 10^^ ergss"^. The redshifts derived from 

the stellar continuum and from the emission lines are 
consistent within the uncertainties. 

Comparison of the spectral continuum emission and 
the rest-frame UV-NIR photometry with SPS models 
implies a stellar mass of ~ 2 X IO^Mq, a reddening 
of Av~ 0.0 — 0.3 mag, a star formation timescale of 
T = 0.2 - 0.3 Gyr, and an age of 1.3 - 2.2 Gyr. The 
results are slightly different for the different SPS mod- 
els and different assumed metallicities. We find a low 
SFR of about 1 — 3 Mq yr~^, implying that the star for- 
mation is strongly quenched and reduced by more than 
99% since its major star formation epoch. If this galaxy 
is typical for quiescent galaxies, quenching of star forma- 
tion is extremely efficient. The constraints on the SFR 
are very tight, and compared to the previously pub lished 
shallow spectrum of this galaxy (jKriek et al.ll2008ah they 
are more stringent by a factor of ~8. 

We do detect a faint Ha emission line, which, based 
on emission-line diagnostics, is not caused by stellar ion- 
ization. However, even if we assume that the Ha line 
was due to star formation, it would result in a low SFR 
of ~ 2 — 4 Moyr~^. One possibility is that obscured 
star formation may have been missed. However, at low 
redshift there is a strong correlation between the dust- 
corre cted luminosity of Ha a n d the bolomet r ic lum inosity 
(e.g. JMoustakas et al. 1 120061 ). iReddv et all ()2006f ) found 
a similar relation at high redshift. Thus, based on the 
faint Ha line, the nature of the line emission, and the 
lack of a detection at 24 /im (I. Cury, private communi- 
cation), we do not expect this galaxy to host obscured 
ongoing starburst regions. 

Although we can confirm our earlier conclusions about 
th e quiescent nature o f the stellar population in 1255- 
(jKriek et al.l l2006bf) . we previously underestimated 
the age. The shallower spectrum showed an apparent 
Balmer break and our fit preferred a younger dusty post- 
starburst galaxy. However, our deeper spectrum shows 
that the galaxy is 1.3 — 2.2 Gyr old, and the optical 
break is dominated by the 4000 A break rather than the 
Balmer jump. A post-starburst (0.5 Gyr) population can 
only account for < 10% of the stellar mass. 

The underestimated age in the shallow spectrum is re- 
lat ed to the overestim ation of the continuum redshift. 
In iKriek et all (|2008aD we found a continuum redshift 
of 2.311qq7, which is within 2a consistent with the 



10 



Kriek et al. 



emission-line redshift of 1255-0. This may imply that 
the uncertainties in our previous work are underesti- 
mated. However, comparison of emission-hne and contin- 
uum r edshifts for 19 emission- line galaxies in lKriek et al.l 
([2008a) demonstrated the reliability of the uncertainties 
of our continuum redshifts, and thus 1255-0 is probably 
"the" largest {2 a) outher. 

The fact that we previously underestimated the age of 
the galaxy shows how difficult it is to estimate ages based 
on "shallow" spectra, let alone broadband photometry. 
Thus, we cannot exclude that the ages of more quiescent 
ga laxies may have be en systematically underestimated 
in lKriek et al.l (|2008af ). The best-fit models imply a for- 
mation redshift of z = 4 — 7 for 1255-0, and this may 
indicate that massive galaxies with strongly suppressed 
star formation exist at even earlier times. 

The rest-frame optical emission-line diagnostics in- 
dicate that 1255 -0 most likely hosts a LINER. In 
iKriek et al.l (pOOl we found that at least 4 out of the 
11 emission- line galaxies in our massive galaxy sample 
at 2.0 < z < 2.7 host an AGN. This study was based 
on relatively shallow spectroscopy (~ 3 hr per galaxy). 
For 1255-0, we did not detect any emission lines in the 
shallow NIR spectrum ()Kriek et al.l l2006br ). Thus, we 
may find actively accreting black holes in more massive 
galaxies, when we obtain deeper data. This suggests the 
possibility that low luminosity AGNs may be very com- 
mon in these objects. 

Although an AGN dominates the line emission, its con- 
tribution to the continuum emission is very low. In the 
H band, in which we measu red the compact size of 1255- 
(jvan Dokkum et al.l l2008t) the contribution is < 8%. 
Thus, an AGN could not be the dominant cause for the 
compact size of this galaxy. A central post-starburst 
population may have a larger effect on our size measure- 
ments. If the stellar population in the center has a lower 
M/L than the outskirts, the size will appear smaller. 
Similarly, if the galaxy experienced recent star formation 
in the outer parts, the M/L will be higher in the center, 
and thus the size will appear larger. Although a post- 



starburst population can only account for < 10% of the 
stellar mass, it can make up ~ 40% of the if-band flux. 
Thus, depending on where this post-starburst popula- 
tion is situated, the light distribution in these compact 
galaxies may be more or less concentrated than the stel- 
lar mass distribution, and our size measurements may be 
over- or underestimated. However, regardless of whether 
or not a post-starburst or AGN is present, this spec- 
trum shows that this compact quiescent galaxy is not 
actively forming new stars, and is primarily composed of 
an evolved stellar population. 

The future generation of NIR spectrographs promises 
higher throughput, in combination with multiplexing. 
This will allow even deeper spectra of larger samples. 
However, the most spectacular results - in particular 
with regard to kinematic measurements of quiescent 
galaxies - are expected to come from NIRSPEC on the 
James Webb Space Telescope (JWST). Without the hin- 
drance of the Earth's atmosphere, it will be possible to 
obtain significantly deeper spectra in only a small frac- 
tion of the exposure time needed for the spectrum pre- 
sented in this paper. 
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APPENDIX 

FITTING AND ASSESSMENT OF SYNTHETIC TEMPLATES 

We developed a custom IDL code named FAST, to fit SPS models to broadband photometry, spectra, or both. 
FAST is compatible with the photometric redshift code EAZY (|Brammer et al.ll2008l ). such that the format of the 
input photometric catalog and filter files is similar. Optionally, the photometric redshifts as derived by EAZY can be 
read in and used by FAST. 

Summarized, FAS T reads in a parameter file, which defines the photometric and spectroscopic catalog s, the SPS 
models (choice from iBruzual & Chariot 2003; Maraston 2005), the IMF(choice from Salpetcr 1955; Kr oupal [200lt 
IChabrienl2003f ). the reddening law, and the fitted grid of stellar population properties. For all properties (age, star 
formation timescale r, dust content Av , metallicity and redshift) the minimum and maximum value, and the step size 
can be defined. Optionally, FAST can calibrate spectra using the broadband fluxes. 

FAST generates a six-dimensional cube of model fluxes for the full stellar population grid and all filters, spectral 
elements, or both. To determine the best-fit parameters, it simply determines the of every point of the model cube. 
In case spectroscopic or photometric redshifts are provided, the redshift will be fixed to the closest value in the grid. 

The confidence levels are calibrated using Monte Carlo simulations. The number of simulations can be defined in 
the parameter file. The observed fluxes are modified according to their p hotometric and spectr oscopic errors in each 
simulation. Optionally, a rest-frame template error function (see EAZY, iBrammer et al.l [20081 ) can be added to the 
broadband photometric fluxes, to account for uncertainties in the models. In case the automatic scale option is used, 
the spectrum is scaled individually for each simulation. By doing this, we incorporate the error on the scaling factor*. 
We determine the best solution for all simulations and define the Xio- (xio- or Xso-) level as the value in the originally 

* Note that we did not apply this method in IKriek et al.l l|2008al '): instead, we arbitrarily increased the uncertainties by quadratically 
adding 10% of the average flux of the spectrum to account for systematic effects. 
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grid that encloses 68% (95% or 99%) of the simulations. The uncertainties on the stellar population properties are the 
minimum and maximum values that are allowed within the Xia (xio- or xia-) threshold. 

In case photometric redshifts (as provided by EAZY) are assumed, the calculation of the confidence intervals is 
slightly more complicated. Both codes use different template sets and thus produce different probability distribution 
functions of z. Unfortunately, there is no perfect method to determine the confidence levels for this case. We use the 
following method, which incorporates the confidence levels of EAZY, and is fast and user friendly. We run all Monte 
Carlo simulations at the best-fit Zphot as derived by EAZY. We determine A^^; the difference between the value 
that encloses 68% (or 95% or 99%) of the simulations and the minimum value at this 2:phot- Next, we return to the 
full gi'id that resulted from FAST. The confidence intervals for a given parameter are the minimum and maximum 
values in the full model cube with a x^ value less than Xmin + '^"■^^ a redshift within the 68% confidence interval as 
provided by EAZY. In essence, we use the EAZY output to limit the redshift range of the solutions allowed in FAST. 
Note that Xmin of ^^e full grid is likely lower than the minimum value of at Zphot- We tested this method using 
a large number of Monte Carlo simulations that are input in both EAZY and FAST (retaining the combinations of 
2phot and perturbed fiuxes). In this case we do not use the simulations for calibration, but derive the uncertainties 
directly from the output. We find good correlations between the uncertainties of the two methods and we believe that 
the FAST confidence levels are robust and reliable. 
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